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Abstract. We present the results of new 3D TIGER 
spectroscopic observations and archived HST/WFPC2 
and NICMOS images of the central region of M 104. The 
[Nll]+Ha images reveal the presence of a nuclear spiral 
structure, and the gaseous kinematics in the central arc- 
second shows evidence for kinematical decoupling of the 
central peak. A straight nuclear dust lane, with a weak 
symmetric counterpart, is seen in the V — I and V — H 
colour maps. These results hint for the presence of a strong 
nuclear bar, that would be located inside the inner Linblad 
resonance of the large-scale bar discussed by Emsellem 

dm. 
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1. Introduction 

This is the third of a series of papers dealing with 
optical observations of the Sombrero galaxy (M 104, 
NGC 4594) . We report here the results of 3D spectroscopy 
of the nuclear region of M 104 with the TIGER instru- 
ment, in the 6750/460 A spectral domain, which includes 
the [Nn]AA6548,6583, Ha and [S n]AA6717,6731 emission 
lines. This 3D dataset is used in combination with archived 
images acquired with the Wide Field and Planetary Cam- 
era 2 (WFPC2) and the Near Infrared Camera and Multi- 
Object Spectrometer (NICMOS) on board the Hubble 
Space Telescope (HST), to study the nuclear regions of 
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this galaxy. This paper is organised as follows: Sect, 
presents the observations and the data reduction; Sect. 
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shows the results; and Sect. [| contains a brief discussion 
of the nature of the observed nuclear structures. 

Throughout this paper, we will u se a distance to M 104 
of 8.8 Mpc (Ciardullo et al. 1993 ), yielding an intrin- 
sic scale of ~43 pc.arcsec -1 . We will also use a systemic 
(heliocentric) velocity V s of 1080 km s _1 , infered from 
the two-dimensional stellar velocity field (Emsellem et al. 
|l996j hereafter Paper 2), and a value of 84° (±2°) for the 
inclination of the galaxy (Paper 2). 



2. Observations and data reduction 

2.1. HST/WFPC2 and NICMOS imaging 

From the Space Telescope/European Coordinating Facil- 
ity (ST/ECF) archive at the European Southern Observa- 
tory, we have retrieved WFPC2 and NICMOS images of 
M 104 in four bands: F547M, F658N, F814W, and F160W. 
Their total integration times were 1340, 1120, 1600, and 
128 seconds, respectively. The data reduction was per- 
formed using the ST/ECF pipeline, with the most recent 
calibration data. 

The images were flux calibrated, rotated to the car- 
dinal orientation (north up, east to the left), and cor- 
rected for geometric distorsion. An emission-line image 
([Nii]+Ha) has been constructed by subtracting the 
F814W image from the on-band F658N image. Colour 
maps, similar to V — I and V — H, were also constructed, 
using the F547M, F814W and F160W images]]. 

2.2. TIGER spectrography 

We obtained a total of 1.5 hour integration on the cen- 
tre of M 10 4 using the TIGER 3D spectrograph (Bacon 
et al. 1995| ) at the Canada France Hawaii telescope, in 
April 1996. The spectral domain (6750/460 A) includes 
the [Nn]AA6548,6583, Ha and [S n]AA6717,6731 emission 
lines. The spectral sampling was 1.5 A per pixel, with a 



The V — I map has alread y been presented by Pogge et al 



( 1999 ) and K ormendy et al. (1996). Note that in Kormendy 
et al. ( 1996] ), the FOS apertures in their Fig. 1 are on the 
wrong side (east and west are reversed). 
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Fig. 1. Top row: [Nn]A6583 line flux (left), heliocentric centroid velocity V (middle, minus V s — 1080 km s" 1 ) and 
dispersion a (right, corrected from the instrumental dispersion) maps, as derived from a single component Gaussian 
fitting of the emission lines in each spectrum of the original (i.e. not deconvolved) TIGER data cube. Black or white 



contours correspond to selected isophotes (2.5, 5, 15, and 45 x 10 18 W m 2 arcsec 
map. Spectra with a [Nn]A6583 line peak intensity < 0.9 x 10~ 19 W m~ 2 A" 1 



2 ) of the [Nn]A6583 line flux 



arcsec 



(i.e. ~ 3(7]] 



5 ) have 



been discarded in the V and a maps. Bottom row: Same as top row, but derived from the deconvolved data cube 



(Richardson-Lucy deconvolution, 40 iterations, see Sect. p.2[ ). Spectra with a fitted [Nn]A6583 line peak intensity < 
" 19 W m~ 2 A -1 arcsec" 2 (i.e. ~ 3<7 n 



0.24 x 10" 



,) have been discarded in the V and a maps. North is up, East left. 



spectral resolution of ~3.5 A FWHM. The data red uctio n 
was performed using a dedicated software (Rousset 1992|) . 
The spatial sampling was 0'.'39 per lens, with a (seeing lim- 
ited) spatial resolution of 0'.'95. Subtraction of the stellar 
continuum was achieved via a library of stellar and galac- 
tic spectra using a procedure which will be detailed in a 
forthcoming paper (Emsellem et al. in preparation). 

To improve the spatial resolution of this data set, we 
have used a Richardson-Lucy algorithm (Richardson 1972; 
Lucy 1974) to deconvolve each velocity slice of our con- 
tinuum subtracted data cube. The point-spread function 
used in this deconvolution was obtained by comparing the 
emission-line image derived from the TIGER data cube 
and the HST emission-line image. We limited the number 
of iterations to 40 and the spatial resolution after decon- 
volution is ~ 0'.'5, with a new spatial sampling of 0'.'2 per 
pixel. 

The emission lines present in our spectra h ave b een 
fitted using the FIT/SPEC software (Rousset |l992[) to 



reconstruct maps of the ionised gas distribution and 
kinematics. Except when explicitely mentioned, we only 
present results obtained by fitting a single Gaussian profile 
for each individual emission line, with all lines constrained 
to share the same velocity shift and width. Excellent fits 
were obtained, even for the nuclear emission-line profiles, 
despite their extended wings (e.g. see Fig. 3 in Kormendy 
et al. |1996|, hereafter K+96). 



3. Results 

3.1. Ionised gas distribution and kinematics 

In Fig. |l|, we present maps of the [Nn]A6583 line inten- 
sity, centroid velocity and velocity dispersion (the latters 
being common to all fitted lines). The centre of the field 
has been defined as the location of the maximum line in- 
tensity, which coincides for all observed lines. The ionised 
gas distribution derived from our deconvolved data cube 
agrees with that obtained in the HST [NnJ+Ha; image 



x [arcsec] 



Fig. 2. Plot of the (heliocentric) ionised gas velocity (mi- 
nus V s = 1080 km s _1 ) along the major axis of the galaxy 
(east is left). Velocities derived from the spectra within 
0'.'2 of the major axis (original dataset) are represented 
by circles, with error bars overimposed. The upper-panel 
shows the location of these spectra on the sky. We also 
display the interpolated velocity profiles, as derived from 
the original (dotted line) and deconvolved (solid line) data 
cubes. 



(Pogge et al. 1999), and exhibits three prominent features: 



a bright core, two roughly symmetric spiral-like structures, 
and a curved extension northwest of the nucleus. 

The spiral structures harbor velocities significantly 
lower than the circular velocities predicted by the multi- 
Gaussian expansion model of M 104 (Paper 2). Similar 
behavior has been noticed fo r gas between 10 and 50" 
from the nucleus (Rubin et al. 1985 ), as well as in nuclear 
FOS spectra (K+96). A cut of the velocity map along the 
major axis of the galaxy (Fig. ^), shows the presence of 
a strong velocity gradient close to the nucleus, with ex- 
trema of +130 and -200 km s -1 , located 0'.'4 (18 pc) east 
and l'.'l (48 pc) west of the centre, respectively. These ex- 
trema are followed by an abrupt decrease of the velocity 
modulus, which reaches minima l'/4 (60 pc) east and 2'.'3 
(100 pc) west of the nucleus, before increasing again. The 
kinematics of the ionised gas inside ~1" (43 pc) is thus 
decoupled from that of the gas in the spiral arms. This 
decoupling is also present in the datacube before decon- 
volution. After deconvolution, the velocity gradient in the 
central arcsecond is larger on the eastern side than on the 
western side. 

The velocity dispersion is centrally peaked, reaching a 
value of 380 km s _1 in the deconvolved data. From their 
high spatial resolution FOS data, Nicholson et al. (1998) 



obtained larger central velocity dispersions, of 540 and 
390 km s _1 for [Nil] and Ha, respectively. However, this 
discrepancy can easily be accounted for by the difference in 
spatial resolution and the fact that we did not attempt to 
model the [N n] and Ha lines separately. Note also that the 
central velocity dispersion for the gas in our deconvolved 
data is similar to the stellar one from K+96 (FOS data). 
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Fig. 3. Top panel: V — I colour map, with selected 
isophotes (5, 20 and 100 x 10~ 18 W m -2 arcsec -2 ) of 
the HST [Nll]+Ha image overimposed. Bottom panel: 
V — H colour map. 



3.2. Colour map and dust distribution 

The colour maps (Fig. ||) outline the (patchy) dust struc- 
tures present in the central region of M 104. As already 
pointed out in Paper 1, there is a straight-line dust lane 
southeast of the nucleus. It has a projected length > 1'.'2 
(~ 50 pc). Taking a minimum inclination of 82° for the 
galaxy and assuming that this dust lane lies in the equa- 
torial plane of the galaxy, this yields a deprojected length 
> 370 pc. The associated E(V — I) reaches 0.2 mag, cor- 
responding to an apparent Ay ~ 0.5 mag (for Ry — 3.1). 
Especially in the V — H colour map (Fig. [| bottom panel), 
there is a hint for the existence of a symmetric, though 
much fainter, extinction feature on the northern side of 
the nucleus. 

Last, the nucleus of the galaxy appears extremely blue 
(local V — I difference of 0.3 mag) compared to its sur- 
roundings. It is unlikely that this feature is an artifact^, 
as already mentioned in K+96. 

4. Discussion 



Recently, Regan & Mulchaey ( 1999j ) argued that nuclear 
bars may not be the primary agent for the fueling of 
active galactic nuclei (AGN) and suggested nucle ar spi- 
rals as an alternative mechanism. Martini & Pogge (1999) 
also analysed visible and near-infrared HST images of 



2 We checked that it is still present after harmonisation of the 
spatial resolutions in the F547M, F814W, and F160W images. 
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a sample of 24 Seyfert 2 galaxies and found that 20 of 
these exhibit nuclear spirals (with only 5 clear nuclear 
bars). Physical processes involved in the formation of such 
spirals are not clear yet, although acoustic instabilities 
have often been mentioned as a possible mechanism in 



non self-gravitating nuclear discs (Montenegro et al. 1999, 
Elmcgrccn et al. 1998, and references therein). In M 104 
(classified as a liner), kiloparsec-scale spiral arms are in- 
deed present but a straight dust lane is also found closer to 
the nucleus: it is a specific signature of s trong bars with 
inner Lindblad resonance (Athanassoula 1992 ). If it ac- 
tually traces a nuclear bar, this dust lane should have a 
symmetric counterpart on the other side of the nucleus. 
There is indeed a hint for such a feature (see Sect. 3.2 ) 
but it is rather weak. However, this weakness could be at- 
tributed to the fact that this second dust lane would be 
located on the far (north) side of the galaxy. 

To test this hypothesis, we computed the effect of the 
presence of two symmetric dust filaments using the lumi- 
nosity density model of M 104 of Emsellem ( |1995| , Paper 1 
hereafter). We assumed the filaments to be in the equa- 
torial plane of the galaxy (10° from end-on), and their 
characteristics were set to ensure that the mean extinc- 
tion for the southern filament was consistent with the 
observations. The model predicts that the apparent ex- 
tinction should rapidly decrease northward, with appar- 
ent A v < 0.03 mag or E(V - H) < 0.026, Of! 6 north of 
the nucleus (for Ry = 3.1). These are upper limits since 
dust scattering and dumpiness would tend to significantly 
reduce these values. For these extinction levels, we indeed 
expect the northern filament to be barely detectable in 
our colour maps. 

The observed kinematics would also fit naturally into 
the strong nuclear bar picture. The fact that overall the 
velocities are small compared to the predicted circular ve- 
locities, obviously argue for the presence of strongly non- 
circular motions. The velocity profile in the spiral features 
could be explained by the combination of the streaming 
motions and projection effects as the spiral curves around 
the nucleus and becomes perpendicular to the line-of- 
sight. If the nuclear bar hypothesis is correct, the existence 
of offset straight dust lanes also requires the presence of 
an inner Linblad resonance (ILR), and an extended Xi 
orbit family (e.g. Athanassoula & Bureau 1999 , hereafter 
AthB99). This would explain the observed kinematical de- 
coupling between the gas in the nucleus and in the spi- 
ral arms, as clearly illustrated in the models of AthB99. 
The observed asymmetry in the central velocity gradi- 
ents agrees also qualitatively with the predictions made 
by AthB99, when dust is included (see their Fig. 10 with 
e.g. ip = 22.5°). We should however keep in mind that 
significant differences exist between the gas distribution, 
as observed in the core of M 104, and as idealised in the 
model of AthB99. 

The presence of a large-scale bar has already been sug- 
gested in Papers 1 and 2 (pattern speed il p of ~ 120 km s _1 



kpc -1 for a distance to M104 of 8.8 Mpc). The mass 
model predicted a strong ILR for this primary bar, located 
roughly 20" from the nucleus. The secondary nuclear bar 
discussed in this letter would then be well inside this res- 
onance. From the extension and orientation of the central 
dust lanes, we can roughly estimate a semi- major axis of 
a = 425 pc (~10"), and an orientation of approximately 
10° from end-on^] for this bar. The location of the coro- 
tation can be estimated from the relation ~ 1.2 x a 
(AthB99, Tl is the Lagrangian radius), yielding a value of 
~12". This is right where the transition region b etwee n 
the inner and outer disks occurs (Seifert & Scorza 1996 ). 
More detailed modelling is needed to accurately estimate 
its pattern speed. 

If there are clearly some hints of the presence of a nu- 
clear bar in M 104, as discussed above, it is also clear that 
additional information are needed to confirm or infirm 
its existence. In particular, the main support for the nu- 
clear bar scenario comes from the presence of the straight, 
southern dust filament, which could actually be much fur- 
ther from the nucleus than we assumed (e.g. if it is out- 
side the equatorial plane of the galaxy) . In the same way, 
the kinematical decoupling of the central regions could al- 
ternatively be due to, e.g. a nuclear keplerian disk. The 
answer to this should come soon, from HST/STIS observa- 
tions of M 104. These high spatial resolution data will al- 
low detailed comparison between the observed PVDs and 
the gaseous kinematics predicted by various models (nu- 
clear bar, keplerian disk). Emission lines in the infrared, 
where dust is less problematic, could also be very valuable 
to understand the gas distribution and kinematics in more 
details. 
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